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The energy momentum tensor of the perfect fluid is a simpHfied but successful 
model of wide practical use in astrophysics. In the previous researches, some im- 
portant information such as the self-energy of the interaction of the particles and 
influence of the gravitational potential are usually ignored or introduced phenomeno- 
logically. In this letter, the analysis and calculation show that, the interactive po- 
tentials of the particles can be described by introducing an extra potential W, which 
acts like negative pressure. For the gases in a star, the state functions are strongly 
influenced by the gravitational potential, and only the temperature T is independent 
variable, but the other state functions such as mass-energy density p and pressure 
P can be expressed as simple functions of T. These relations can provide some new 
insights into the dynamical behavior and structure of a star. 
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In astrophysics, the energy momentum tensor of the particles in a star is usually described 
by the Einstein's perfect fluid model 

T^" = {p + P)U^U'' - Pg^". (1) 

Theoretical analysis and experiments all show (1) is a successful approximation. Due to the 
effort of several researchers such as Israel and Stewart [1, 2], Carter [3] and Lichnerowicz[4], 
the energy momentum tensor T'^'^ has been developed into a treasure box for classical fluid 
theory, which contains not only energy and momentum, but also heat flux, spatial stress 
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and viscosity [5, 6]. Furthermore, the relativistic hydrodynamics with variational principle 
and Cartan's exterior algebra was elegantly established [7]. 

However, some important information is ignored in the perfect fluid model (1). One is the 
self-energy of the internal interactions of the particles such as electromagnetic potential and 
nonlinear potential of a spinor, the other is influence of the gravity on the motion of particles. 
In this letter, we discuss the two realistic problems in detail. The discussion shows that the 
self-energy of interaction can be merged into the mass-energy density, and an extra potential 
term W should be introduced into T'^'^ which acts as negative pressure. The meaning of 
mass-energy density p should be redefined to include the contribution of the self-energy. The 
gases in thermal equilibrium inside a star is a state automatically formed, so the equation 
of state can not be arbitrarily given, and only the temperature is the independent variable. 
This property is much helpful to study the interior behave of a star, and leads to an elegant 
and natural model for stellar structure [(S]. Besides, the analysis also shows the consistence 
and connection among thermodynamics, classical mechanics and field theory. 

In conventional thermodynamics, some concepts are usually not explained definitely, e.g. 
even the temperature T was not clearly defined for large energy range. This is inconve- 
nient for theoretical analysis and easily leads to confusion. In the case of the monatomic 
gases, considering that the thermal motion of particles in the comoving coordinate system 
is fundamental, so we define the temperature T of the gases as follows 

N 

3" 3N 



kT=-K = --yKn, (2) 



where k is the Boltzmann constant, K the mean shifting kinetic energy of particles in a 
micro comoving volume element V, N the particle number in V, Kn the kinetic energy of 
n-th particle. (2) is the most popular description in the usual textbooks. 
Assuming the distribution of kinetic energy of the particles is given by 

''-/(^)i^ (3) 

which satisfies 

l*oo poo q roc o 

dV = l, / KdV = -kT, / K^dV = —{kTY, (4) 
Jo Jo 2 Jo 2cr 

where a is a constant related to the concrete distribution function for given particles. In usual 

statistical mechanics, the distribution functions are usually defined in the phase space of 
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coordinate and momentum of the particles, but such selection is sometimes inconvenient for 
calculation of state functions, especially in the case of relativistic gases. Since the following 
discussions have not any explicit relation with distribution function /, and at most use the 
second order moment, so the kinetic energy distribution (3) is much convenient and the 
quantities defined by Eq(4) are enough. The following examples show how a varies with 
distribution /, 

dV = ( ) Y '^^'j^' ^ " 5' (^^^^^^^ distribution), (5) 

40 f K\\'^ dK . 15 



i'+177^"H^JJ ln(217/40)fcT ' "=38' 
/ 2K Y'' 2{n-l)dK n-3 

V 3(r2 - 2)kT ) 3{n-2)kT' ~ 3(n-2)' ^ ^ 

, f K 3\ dK 2 , , . . r ^ 

dV = 2J "fcT ' " = 3' ^ ^^^^^ bound of a). (8) 

For the quantized nonlinear fields such as spinors, in microscopic view, the energy mo- 
mentum tensor can be expressed by the following classical approximation [9, 10, 11] 



T^^ = J](m.« + w.gnS'iS - X,)^l - vl (9) 

k 

where is the proper mass of the /c-th particle, Wk is the proper potential, 4-vector 
velocity, Vk the speed, Xk{t) the central coordinate. In the case of ideal gas, we have 
Wk = 0, (V/c). The physical meaning of Wk can be explained by the following example. 
Consider the following Lagrangian of self-interaction 

Ca = Y^ CA\k, CA\k = ^{d^A.^kd^Al - sM^A„|fc) + e^g^A^^, (10) 

k 

where stands for the self potential of the k-th spinor, s = corresponds to the electro- 
magnetic field. Making classical approximation, we generally have[ll, 12, 13] 



- J2i(<+P>'k''^k-P'>^9n^l-vl5'ix-Xk), (11) 

k 

where Xk is the central coordinate of the particle, m'^ stands for the proper energy of A'j^ 
in mean sense, and p'j. for the internal pressure. By (11), the parameters in the central 
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coordinate system are calculated by 

m'f. = - CA\kdx = l-ek / q^^Au^d^x > 0, (12) 

Jif3 ^ Jffi 

p'k = I (rn', + jjVAlVA^\u + ACA\k)d?x^ = (^m', + s'AtA^\,d'x^ , (13) 

in which we have used doA°' = and —AA'^ + s'^A'^ = eq^ in the central coordinate system. 
From (13) we learn p'f. < 0. So different from the photons which have = |pph; the static 
field A'^ provides negative pressure. By (12) and (13), we find (9) is of general meaning, in 
which rrik can be explained as the total mechanical mass of the k-th particle including the 
self-energy of interaction, p'^ < is the total potential acting as negative pressure, which 
leads to Wk > in (9). 

In order to calculate the statistical expectation value of Eq(9), denote 

P^-hTl + Ty^ + T%), K, = ^^-mk, (14) 

then we have 




In Eq(15) and (16), we can directly solve the mean value of K^, rrik as the function of kT and 
fh by (4), but the rest terms include the factor -^/l — v^, which is difficult for estimation. 
We resolve the problem by the following trick, which naturally includes the influence of the 
gravitational potential. 

To research the thermodynamical properties of gas, we use piston and cylinder to drive 
gas. In astrophysics, we have more ideal piston and cylinder, that is the spacetime with 
Friedmann-Robertson-Walker(FRW) metric, which is absolutely adiabatic and reversible. 
The FRW metric drives the gases homogeneously expanding and contracting as the scale 
factor a varies, and the results have generally meanings according to the principle of equiv- 
alence. 

The line element of the FRW spacetime is given by 

ds^ = de - a'it) [dr^ + h^{r)de^ + h\r) sin^ Odip^) , (17) 
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where the specific form of h{r) is not important. Solving the geodesic we have the following 
relations[10, 14, 15] 



kT = , , r^r-^, , (19) 



^kT{kT + 2am) ' a(a + + fe^) ' 

where 6^, ^ are determined by initial velocity. Since we firstly solve the micro motion of 
particles and then solving the expectation value of a variable, the results should be the same 
as that we firstly do statistical average for the variable and then solve its macro changing 
value, so we have the following calculation. By (18), we have relation 

~ da Id 



da — i;2 rrik da 
Then by (19), we get 



a 



{Kk + rrik)]. (20) 



^ E - / y^^^ + 'n.)^v) ^ r-n - 5^^. (21, 



where the mean parameters are defined by 

N ^ N 



k k 



In the case Wfc > 0, by Eq(22) we have 

E J ^/^kdV > 0, a-^< ^. (24) 

k 

For the same kind particles, by w = firh, we get o" < |- 

Similar to thermodynamical laws have not relation with piston and cylinder, Eq(21) 
and (22) also have not any relation with the 'piston' a(t), so they are independent of any 
dynamical process, and are generally valid for ideal gases in the comoving coordinate system. 

We define the mass-energy density p and pressure P of the particles in comoving coordi- 
nate system as usual [9] 

'^\/l-vt i'-Tx/X-vl 
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Then, in mean sense we get 

P=l^J dV jj£^{Kk + mu)5\x-X,)dV = Q(^ + ^-^^, (26) 

where q = is the static mass density. By Eq(15), (16) and (21), (22), we find the 
expectation value of Eq(9) should take the following form 

T^"" = {p + Pp^'U" + {W - P)g^"', (27) 

where is the mean velocity satisfying U^U^ = 1, and W is an ignored term corresponding 
to potentials Wk, which acts as negative pressure. 

Then according to (21) and (22), Eq(26) and (27) give 

.^i/.p/n--P^.(|-,^), (2s) 

P^l/.p/p.V^ + ^r^.f (l-^^-^). (29) 

In case of ideal gases or spinors, Eq(29) can be also derived according to the principle of 
action[14]. In the ultra-relativistic case kT ^ rh, we have P ^ , which is obviously 
reasonable in physics. 

W can be explained as potentials. If Wk = 0, we get 1^ = 0, and the other results reduce 
to the relations for ideal gas[14, 15]. So the introduction of W provides an entry to deal 
with the self-energy of internal interaction potentials of the particles. 

The equations (26), (28) and (29) are generally valid relations for any ideal gases in local 
equilibrium. By (19), noticing g = and oc oc [kT(kT + 2crm)]~t, we get the 
following state functions in dimensionless form 

P = mWJ + 2<7)]}. (30) 

p = £)(l + |j) = i)„|J(J + 2ff)ll(l + ^J), (31) 

^-^{^-^)-^«W^2.)]i{,-J£^), (33) 

where Qq is a constant depending on parameters (m, b, a, /i, fi), and 
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Surprisingly, under the constraint of gravity, among the state functions only the temperature 
is independent. 

Equations (31)- (33) are based on the assumption that the particles move along geodesic. 
But the analysis in [11] shows that the central trajectory of a spinor with self interactions, 
such as nonlinear potential and electromagnetic field, will slightly depart from the exact 
geodesic. So in strict sense, we should solve the exact trajectories to replace the equations of 
geodesies (18) and (19). This is a complicated problem. However, in the case < w;^ <^ m^, 
the above results are good approximations to high precision. Moreover these results can be 
modified by the following treatments. 

In (31)-(33), we treat density q or equivalently the volume V as an independent state 
function, and then derive the function q{J) according to the continuity equation. In the 
comoving coordinate system with the following Gaussian type metric [9] 

g^y = diag {goo, -gab) , (a, 6) G {1, 2, 3}, (35) 

where gab is the spatial metric, for energy- momentum tensor (27), we check the following 
Gibbs-Duhem's law 

6Q = d[{p + W)V] + {P-W)dV, (36) 

where SQ denotes the heat received by the particles, and 

V = v^AxA|/A2;, (g = det{gab)) (37) 

is the micro spatial volume occupied by given particles. 

We trace the motion of these particles. In the adiabatic process, we have SQ = 0. Then 
by U^T^^'^ = 0, we get the continuity equation for Eq(27) as 

W'd^ip +W) + {p + P)U''.^ = 0. (38) 

Denoting the proper time by dr, then ^ = U^d^, Eq(38) becomes 
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where g = det{g^,y) = goo9- In the comoving system, we have = {^^/g^, 0,0,0) and 
dr = ^/gmdt, then we get 

^d^,U^ + -^^/g^= ^/g^dt-^ + -^dty^=0. (40) 
dr Jg^ Jg^ 



Multiplying Eq(39) by the comoving volume element AxAyA^ for the given particles, 
(39) gives (36) in the case 5Q = 0. 

One may oppugn the generality of Eq(36) due to the collision of particles. As shown 
in [14], Eq(36) holds in mean sense for any process with elastic collisions, which can be 
understood as relay race of particles. (36) include the gravitational potential of the particles, 
which is different from but consistent with the energy- momentum conservation law T'^J^ = 0. 
For example, T^^ = holds for a little solid ball, in which the molecules are all fixed at 
definite position by electromagnetic interaction, and then the influence of the gravity on the 
thermal movement vanishes. However, in astrophysics the gravity is dominant inside a star, 
where the gravity drives the particles moving rapidly, and such motion should be exchanged 
and balanced with thermal motion of the particles. This is the underlying and subtle relation 
between gravity and equation of state, which is ignored before. If the gravitational potential 
is taken into account, we get elegant and singularity- free models for the stellar structure[8]. 

In the case of 6Q ^ 0, (36) provides an entry to deal with the radiation and heat produced 
by nuclear reaction inside a star. Obviously, 6Q takes the form 6Q = {gq{J) — -Eradiation)'^''"- 
Then in the equilibrium process, we also have 6Q = 0. 

Now we derive the relation g{J) from (36) for the equilibrium process with 5Q = 0. For 
clearness, we use the the static mass density g = to replace V. Substituting Eq(26), 
(28) and (29) into Eq(36), we get dimensionless differential equation 

Idg 3[(J + (t)2 -/xa^] 



gdJ {J + a)[.P + BJ -2^ia\ 

and the solution 



(41) 



g = g,(J + a)WTT (j2 ^ _ 2a/i) ^ ^ [ 2J + B + a ) ' ^^^^ 

where parameters are defined by 

B = {2 + 3fi)a - 2fi, A = ^^a^i + B"^, a = 3(1 + 4/i)[2/i - 3^- (43) 

(42) shows how the internal potentials Wk influence the mass density. (42) reduces to (31) 
if = /i = 0. 
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In summary, the above calculations work out all relations in the energy-momentum tensor 
of the monatomic particles, which include the classical approximation for the self-energy of 
interaction and the influence of the gravitational potential. Although the realistic process 
inside a star is complicated, the above treatment includes the dominant parts of these effects. 
Noticing the effectiveness of the assumptions, the results will provide some new insights into 
the interior behavior and structure of the stars. In addition, we list some interesting prop- 
erties of the state functions derived above, which show the profound relationship between 
the empirical laws and the theoretical concepts: 

1. In the case of ideal gas with w = ai low temperature, (31)-(32) gives the equation 
of state for the adiabatic monatomic gas 

P=pJ=P^pl^ (J<(j, or fcT<m), (44) 

which is identical to the empirical law in thermodynamics. When J ^ 1, we have 
P — > |p, so the adiabatic index is not a constant for large range of temperature due to 
the relativistic effect. These results show the validity of (31)- (33) and the consistence 
with normal thermodynamics. 

2. By (31), letting J — > cxd or m — >■ 0, we get the Stefan-Boltzmann's law p oc T^. 
This means that the above results automatically include photons, and the Stefan- 
Boltzmann's law is also valid for the ultra-relativistic particles. 

3. In general relativity, all processes occur automatically, and is independent of any 
practical process, so it is a constant representing the strength of gravity constraining 
free particles. Of course, Qq is related to the property of particles. Furthermore, 
equation of state (31)-(32) provides a self consistent theory for stellar structure in 
thermal equilibrium [8]. 

4. In the case w > 0, the motion of the particles with internal interaction will slightly 
deviate from the geodesic [11], so consequently the above relations have a very tiny 
violation. By (42) we find, g = leads to J = ~ p, which means the zero 
temperature can not reach. This result may be caused by the deviation of motion 
from geodesic Eq(18). 
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